Fhree-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD

Based on the successive iteration in the Taylor series expansion method, a three-point explicit compact difference scheme with arbitrary order of accuracy is derived in this paper. Numerical characteristics of the scheme are studied by the Fourier analysisl Unlike the conventional compact difference schemes which need to solve the equation to obtain the unknown derivatives in each node, the proposed scheme is explicit and can achieve arbitrary order of accuracy in space. Application examples for the convectiondiffusion problem with a sharp front gradient and the typical lid-driven cavity flow are given. It is found that the proposed compact scheme is not only simple to implement and economical to use, but also is effective to simulate the convection-dominated problem and obtain high-order accurate solution in coarse grid systems.     

复杂无粘流场数值模拟的矩形/三角形混合网格技术

建立了一套模拟复杂无粘流场的矩形／三角形混合网格技术，其中三角形仅限于物面附近，发挥非结构网格的几何灵活性，以少量的网格模拟复杂外型；同时在以外的区域采用矩形结构网格，发挥矩形网格计算简单快速的优势，有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足．混合网格系统由修正的四分树方法生成．将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算，消除了全矩形网格模拟曲边界的台阶效应，同时保证了网格间的通量守恒．数值实验表明本方法在模拟复杂无粘流场方面的灵活性和高效性．     

Analysis and Modeling of the Turbulent Diffusion of Turbulent Kinetic Energy in Natural Convection

Buoyant flows often contain regions with unstable and stable thermal stratification from which counter gradient turbulent fluxes are resulting, e.g. fluxes of heat or of any turbulence quantity. Basing on investigations in meteorology an improvement in the standard gradient-diffusion model for turbulent diffusion of turbulent kinetic energy is discussed. The two closure terms of the turbulent diffusion, the velocity-fluctuation triple correlation and the velocity-pressure fluctuation correlation, are investigated based on Direct Numerical Simulation (DNS) data for an internally heated fluid layer and for Rayleigh–Bénard convection. As a result it is decided to extend the standard gradient-diffusion model for the turbulent energy diffusion by modeling its closure terms separately. Coupling of two models leads to an extended RANS model for the turbulent energy diffusion. The involved closure term, the turbulent diffusion of heat flux, is studied based on its transport equation. This results in a buoyancy-extended version of the Daly and Harlow model. The models for all closure terms and for the turbulent energy diffusion are validated with the help of DNS data for internally heated fluid layers with Prandtl number Pr = 7 and for Rayleigh–Bénard convection with Pr = 0.71. It is found that the buoyancy-extended diffusion model which involves also a transport equation for the variance of the vertical velocity fluctuation gives improved turbulent energy diffusion data for the combined case with local stable and unstable stratification and that it allows for the required counter gradient energy flux.     

高精度非线性格式WCNS的分析研究与其应用

首先将Fourier方法推广于高维方向研究了五阶精度WCNS的特性,并与其他高阶格式进行比较。分析结果表明WCNS的高精度特性普遍接近甚至好于迎风偏置五阶显式格式EUW-5与Pade′标准格式。然后开展了WCNS的应用研究,采用高效率的WCNS-E-5数值模拟了含强激波的高维复杂流场。算例包括二维高超声速边界层对自由流扰动波的吸收问题以及三维球头绕流问题。计算结果反映出WCNS-E-5对激波等间断的良好捕捉能力,图像清晰光滑,数据准确可靠。     

An explicit algebraic Reynolds stress and heat flux model for incompressible turbulence: Part I Non-isothermal flow

Tensor representation theory is used to derive an explicit algebraic model that consists of an explicit algebraic stress model (EASM) and an explicit algebraic heat flux model (EAHFM) for two-dimensional (2-D) incompressible non-isothermal turbulent flows. The representation methodology used for the heat flux vector is adapted from that used for the polynomial representation of the Reynolds stress anisotropy tensor. Since the methodology is based on the formation of invariants from either vector or tensor basis sets, it is possible to derive explicit polynomial vector expansions for the heat flux vector. The resulting EAHFM is necessarily coupled with the turbulent velocity field through an EASM for the Reynolds stress anisotropy. An EASM has previously been derived by Jongen and Gatski [10]. Therefore, it is used in conjunction with the derived EAHFM to form the explicit algebraic model for incompressible 2-D flows. This explicit algebraic model is analyzed and compared with previous formulations including its ability to approximate the commonly accepted value for the turbulent Prandtl number. The effect of pressure-scrambling vector model calibration on predictive performance is also assessed. Finally, the explicit algebraic model is validated against a 2-D homogeneous shear flow with a variety of thermal gradients. Dedicated to the memory of the late Professor Charles G. Speziale of Boston University     

Supersonic Inviscid Corner Flows with Regular and Irregular Shock Interaction

The results of a numerical and analytical investigation of steady-state supersonic inviscid flows in corners formed by intersecting compression wedges are presented. The flows considered are symmetric about the corner bisector. The distinctive features of flow pattern formation related with the reflection of wedge-generated shocks from the bisector plane are studied. The wedge angles at which transition from regular to irregular shock reflection occurs are determined both numerically and analytically using the criteria available for plane flows; the data thus obtained are found to be in agreement. Flow patterns with irregular shock reflection, namely, single, transitional, and double Mach, as well as von Neumann reflection, are identified; they are similar to the known types of reflection for plane quasi-steady-state flows. Varieties of these types not observed in the plane flows are found to exist. The effects of the angle of inclination of the plane surfaces of the corner to the freestream direction, the sweep angle of the leading edges, and the dihedral angle are investigated. Some previously unknown parameters of corner configurations for which transition may occur in accordance with the von Neumann criterion are determined.     

Some remarks about the resolution of high velocity flows near low densities

High velocity flows which are exposed to strong rarefaction waves and creating low densities regions in it present difficulties and inaccuracies for numerical resolution. In particular, variables related to the internal energy are wrongly evaluated. Use of classical schemes solving the Euler equations in conservative variables introduces significant errors in the determination of temperature. We recommend to employ a non-conservative formulation of the energy equation. Results found to be more accurate in using the present internal energy formulation. In order to have the formulation available for both shock and strong rarefaction waves, we propose a hybrid formulation of conservative and non-conservative ones, depending on a shock indicator. The results are compared with exact solutions and show a significant improvement of the accuracy. The method is then extended to two-dimensional cases. Received 28 March 1997 / Accepted 18 June 1997     

Droplet break-up through an oblique shock wave

The interaction of a two-phase flow with a wedge where a stationary shock wave is initially settled is studied in a two-dimensional configuration. Before the introduction of the dispersed phase, the flow around the wedge is a supersonic one phase flow such as an attached stationary shock wave is present. Then, the dispersed phase is introduced upstream the initial position of the stationary shock wave. The purpose of this study is to point out two-phase and droplets break-up effects on the oblique shock wave. The two-dimensional equations are solved by a TVD scheme where fluxes are computed by using Riemann solver for the gas phase equations and also for the dispersed phase equations wich is an original approach due to the authors (Saurel et al. 1994). In addition to drag forces and heat and mass transfers, the process of droplets fragmentation based on the particle oscillation is considered. Accepted April 28, 1995     

群桩作用机理的数值试验研究

 在分析桩土支撑体系及其相互作用关系的基础上，利用有限元分析软件ANSYS建立了群桩 体系计算模型. 通过对均质土体例题计算结果分析，可以看出利用该建模方法与传统解答有 很好的一致性. 根据辽宁工程技术大学实验馆场地土实测资料，利用该建模方法，求解出了承台与桩的 荷载分担比，对安全、经济地进行群桩基础设计具有重要指导作用和应用价值.     

稀薄流到连续流的气体运动论模型方程算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数可描述不同流域气体流动特性的气体分子速度分布函数的简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。推广应用改进的Gauss-Hermite无穷积分法和华罗庚－王元提出的以单和逼近重积分的黄金分割数论积分方法等，对离散速度空间进行宏观取矩获取物理空间各点的气体流动参数，由此发展一套从稀薄流到连续流各流域统一的气体运动论数值算法。通过对不同Knudsen数下一维激波管问题、二维圆柱绕流和三维球体绕流的初步数值实验表明文中发展的数值算法是可行的。